% 创建时间：2025-06-20 22:57:24
% 作者：兪馳
clc;
clear;

A=[1 2 3 4;
    4 3 2 1;
    1 3 2 4;
    4 1 3 2]; % 定义矩阵 A
B=[5 1;
    4 2;
    3 3;
    2 4]; % 定义矩阵 B

x=inv(A)*B, % 使用数值计算求解线性方程组 Ax = B，得到解 x
e1=norm(A*x-B), % 计算数值解的误差
x1=inv(sym(A))*B, % 使用符号计算求解线性方程组 Ax = B，得到精确解 x1
e2=norm(A*x1-B) % 计算符号解的误差